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PREFACE 



This manual contains the user's instructions for the FORTRAN compiler 
associated with the Heterogeneous Element Processor (HEP) . It is directed 
to the FORTRAN programmer and the programming staff responsible for main- 
taining the compiler. Its purpose is to explain the HEP version of the 
FORTRAN language, the various inputs and outputs of the compiler, and to 
give some examples. 



SECTION I - FORTRAN SOURCE PROGRAM FORM 



1. Introduction 



A FORTRAN source program consists of one main program and any number of 
subprograms. The main program and subprograms are made up of statements 
using the FORTRAN character set. 

1.1 Character Set 

The character set has two subsets: alphanumeric characters and special 
characters. 

ALPHANUMERIC CHARACTERS: 

Letters (A-Z) and 
Digits (0-9) 

SPECIAL CHARACTERS: 

Blank 
- Equals 
+ Plus 

Minus 
* Asterisk 
/ Slash 

( Left parenthesis 
) Right parenthesis 
, Comma 

Decimal point 
$ Dollar sign 

' Apostrophe (single quote) 

Blanks may appear anywhere in a source program. They are significant only 
in Hollerith constants and format specifications. 
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1.2 Source Statements 

Standard FORTRAN statements are accepted as formatted lines (or records) 
of eighty or fewer characters. Each line is divided into four fields: 



COLUMNS 


FIELDS 


1-5 


Statement label 


6 


Continuation indicator 


7-72 


Statement 


73-80 


Identification 



The identification field is ignored by the FORTRAN compiler and is provided 
for the convenience of the programmer. 

1.2.1 Statement Labels 

The statement label is made up of digits placed anywhere in columns 1-5 of 
the initial line of a statement. Blanks and leading zeros are ignored. 

1.2.2 Statements 

A statement consists of an initial line and any number of continuation 
lines. 

An initial line is a line that is neither a comment line nor an END line 
and contains either a blank or a zero in column 6. A continuation line 
is not a comment line and contains any character other than blank or zero 
in column 6. An END line is a line containing an END statement, which 
cannot be continued. 
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EXAMPLE: 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 ... 






2 








X 
X 


A 
C 

+ 


D 


B 


•f 







The first line of this example is an initial line.^ The second and 
third lines are continuation lines. The statement label is 200. 
The statement is equivalent to: 

A=B+C+D 

1.2.3 Comments 

Comment lines have the character C or an asterisk (*) in column 1. Comments 
are for the convenience of the programmer and permit program descriptions 
to be embedded in the code. Comments do not influence the compiler pro- 
gram except that they are printed in the program listing. Columns 2-72 
may be used in any desired format. 

A comment line can be followed only by another comment line or the initial 
line of a statement. 

1.2.4 Blank Lines 

Blank lines may be included in the source text. Columns 1-72 must be blank. 
A blank line can be followed only by a comment line, another blank line or 
the initial line of a statement. 

1.2.5 Source Statement Ordering 



Table 1.1 shows the sequence in which the source statement groups of each 
program must be written. Within each group, the statements may be written 
in any sequence. 
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DATA and CONST statements may appear anywhere after Group 2 and before 
Group 6, but must appear after any declarations (COMMON, DIMENSION, or 
type) affecting the variables to be initialized. 

FORMAT statements may appear anywhere before Group 6. 
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Table 1.1 



SOURCE STATEMENT GROUP 

PROGRAM 1 

FUNCTION 

SUBROUTINE 

BLOCK DATA 

IMPLICIT 2 - 

COMMON 3 

COMPLEX 

DIMENSION] 

DOUBLE [PRECISION] 

EQUIV[ALENCE] 

EXTERNAL 

INTEGER 

LOGICAL 

REAL 

Statement Function 4 

Assignment 5 

ASSIGN 

Assign and Unconditional GO TO 

BACKSPACE 

CALL 

Computed GO TO 

CONTINUE 

CREATE 

DO 

ENDFILE 

IF 

PAUSE 

PURGE 

READ 

RESUME 

RETURN 

REWIND 

STOP 

WRITE 



END 
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1.3 D in Column 1, Conditional Compilation 

All source lines with the character D in column one are compiled as ordinary 
source statements when the compiler is run with the Conditional Compilation 
option selected. When Conditional Compilation is not selected, these lines 
are treated as comments. 

1.4 COPY Statement 

The copy feature enables parts of the program to be stored in more than 
one file. This statement appears as follows: 

COPY name 

where "name" is a file-path name. The name must be completely contained 
on one line (not extended across continuation lines) and must not contain 
blanks . 

A program may contain any number of COPY statements, but they must not 
be nested deeper than three. 

The contents of the named file are inserted into the source program so that 
the first record of the file is the next line after the COPY statement. 
Thus COPY statements may be labeled and referenced the same as CONTINUE 
statements. A COPY statement must not precede Group 1 statements, and 
must not follow a Group 6 statement. 
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SECTION II - DATA 



2. Introduction 



Data is represented as constants and variables. A constant is a quantity 
whose value is explicitly stated; a variable is a quantity whose value 
may change. Each variable is referenced by an identifier that symbolically 
identifies the variable. 

2.1 Identifiers 

Identifiers are used to give names to: 

- Variables 

- Subprograms 

- Common blocks 

An identifier is a string of alphanumeric characters, the first of which 
must be alphabetic. Any number of characters is allowed in an identifier 
but only the first eight are recognized, Certain identifiers may be 
preceded by a dollar sign ($) (see section 2.2). The dollar sign is not 
included as one of the significant characters of an identifier. 

EXAMPLES : 

X15 

PERMUTATION 

STRAIN 

2.2 Variables 

Variable names are identifiers that represent quantities which may assume 
a number of different values. Each variable has an associated data type: 
integer, logical, real, double precision or complex. The variable may 
only assume values of that data type. 



-7. 



A variable may be either synchronous (the FORTRAN standard) or asynchronous 
An asynchronous variable is actually a pair: a standard variable and an 
access control. Assigning a value to the variable causes a wait until 
its access state is EMPTY and sets its access state to FULL with the new 
value. Conversely each use of any asynchronous variable causes a wait 
until the access state of the variable is FULL, and leaves the access 
state EMPTY. The access state may be interrogated with an intrinsic 
function (FULL or EMPTY) to determine its state. An asynchronous variable 
is indicated by a dollar sign ($) preceding the first character of the 
identifier. Also, the same variable must not be used as both a syn- 
chronous and an asynchronous variable; i.e., $A and A must not both 
appear in the same program unit. 

Variables may be simple or arrays. 

2.2.1 Simple Variables 

A simple variable name identifies the location in which a single variable 
value can be stored. 

EXAMPLES: 

N 

X4 

$VALUE 
LOOK UP 

2.2.2 Arrays 

An array variable name identifies an ordered set of data having one, two, 
or three dimensions. Each element of an array is referenced by the array 
name followed by a set of subscripts. 

An array has an associated data type the same as a simple variable. An 
array must be declared by an array declarator which establishes the number 
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of dimensions and the size of each. The data type of an array is identified 
the same as for a simple variable. 

The declaration of an array type may occur separately from the declaration 
of its dimensions. 

2.2.2.1 Array Storage Allocation 

The multi-dimensional arrays declared by the programmer are assigned to the 
one-dimensional computer memory in such a way that the left-most subscript 
varies most rapidly and the right-most subscript least rapidly. 

EXAMPLE: 

Allocation of array K of three dimensions of two elements each. 



Memory Sequence 


Array Element 


1 


K(l,l,l) 


2 


K(2,l,l) 


3 


K(l,2,l) 


4 


K(2,2,l) 


5 


K(l,l,2) 


6 


K(2,l,2) 


7 


K(l,2,2) 


8 


K(2,2,2) 



2.2.2.2 Array References 

Most FORTRAN statements operate on only one element of an array at a 
time. A member of an array is referenced in the form: 

array(Si,S 2 ,S 3 ) 



where "array" is the name of the array, and the S, are suscript expressions, 
A subscript expression may be any expression (see section 4.1). The expres- 
sion is converted to type integer after evaluation. 
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EXAMPLES : 

Y (1) 

STATION (K) 

Q (LINE (N,X) + RHO,N) 

The value of a subscript must be within the limits specified for the array, 
The number of subscripts must equal the number of dimensions specified 
for the array. 

2.3 Data Types 

The following data types are defined: 

Integer 

Logical 

Real 

Double precision 

Complex 

The name assigned to a variable is associated with a data type either 
implicitly or explicitly. 

The explicit type declaration statements (section 3.3.1) assign the type 
explicitly. 

If the data type of a variable is not declared explicitly, the compiler 
types it implicitly according to the following conventions. 

- If the data identifier begins with one of the letters I, J, K, 
L, M, or N, the type is INTEGER. 

- Any other first letter implies REAL. 

The IMPLICIT statement may be used to alter the above conventions and/or 
introduce similar rules for data types other than integer and real. See 
section 3.3.2. 
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For an asynchronous identifier, the first character after the dollar sign 
is used for implicit typing. 

2.4 Constants 

The data type of a constant is determined by its form. 

2.4.1 Integer Constants 

An integer constant is a signed or unsigned string of decimal digits. It 
consists of up to 19 decimal digits in the range -2 63 Cw9.2xl0 18 ) to 
+2 G3 -l(^9.2xl0 18 ). 

The internal representation of an integer is a full word (64-bit) two's 
complement number. However, integer division produces a result accurate 
to only 56 bits. 

EXAMPLES: 

-1 

1234567890 



2.4.2 Real Constants 

A real constant is a signed or unsigned string of decimal digits that in- 
cludes a decimal point and/or an exponent. Any number of digits may be 
included, but only the first fifteen are significant. A real constant has 
one of the following forms; 
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FORMS : 

±i. 

±i.E+e 

±.i 

± . iE+e 

±i.i 

±iE+e 

±i.iE+e 

where i is a string of digits and e is a 1-digit or 2-digit exponent to the 
base 10. The plus (+) character is optional. The magnitude of a real con- 
stant is in the approximate range of 5.4x10 79 to 7.2xl0 75 . 

EXAMPLES : 

3.1415 

+0.031415E+2 
.031415E2 
31.415E-1 



Internally, real data is stored as follows: 



1 




78 




63 


s 


c 




M 





where S is the sign of the mantissa (+ = 0, - « 1) . 
C is the base sixteen characteristic plus 64. 
M is the mantissa (56 bits) , represented as a sign-magnitude 
number . 



The mantissa is hex-normalized. 
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2.4.3 Double Precision Constants 

A double precision constant is a signed or unsigned string of decimal digits 
that includes an optional decimal point and an exponent. Any number of 
digits may be included, but only the first twenty-nine are significant. A 
double precision constant has one of the following forms: 



FORMS: 



±i.iD±e 
±i.D±e 
± . iD±e 
±iD±e 



where i is a string of digits and e is a 1-digit or 2-digit exponent to 
the base 10. The plus (+) character is optional. The magnitude of a 



79. 



double precision constant is in the approximate range of 5.4x10 to 



7.2x10 



75 



EXAMPLES : 

.31415927D1 
3.1415D0 
31415. 93D-04 

Internally, a double precision constant appears as a real constant followed 
by an integer extension of the mantissa: 



Si 


Ci 


Mi 



M 2 
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S x - Sign 

Ci - Characteristic 

Mi - Mantissai (56 bits) 

Mi - Mantissa2 (64 bits) 

2.4.4 Complex Constants 

A complex constant appears as an ordered pair of real constants enclosed in 
parentheses: 

(ri, r 2 ) 
The real part is ri and the imaginary part is r 2 . 

EXAMPLES: 

(0., -1.) 
(5.2, 2.6) 
(-3.E12, .017E-16) 

Internally, complex data is stored as two real constants (see 2.4.2). 

2.4.5 Logical Constants 

A logical constant is a truth value. 

.TRUE, or .T. .FALSE, or ,F. 

A logical constant occupies one word which is either zero (.FALSE.) or non- 
zero (.TRUE.). 

2.4.5 Hollerith Constants 

A Hollerith constant is a string of characters which is represented in- 
ternally by 8-bit ASCII codes. A Hollerith constant may appear wherever 
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an expression may appear, however it must not appear within an expression. 
It is treated as an integer constant. There are four forms of Hollerith 
constants: 



FORMS: 

nHCiC2 . . »C 



n 



*CiC2 . . .C 



n 



nRCiC2 . . . C 
nLCiC2 • . »C 



n 



n 



where n is an unsigned integer constant and C. are characters. 

Internally, the characters are packed eight per word. The two forms on 
the left are equivalent and imply left justification with trailing blanks 
(if necessary to pad n to a multiple of 8) . The R and L forms indicate 
right and left justification respectively, however these two forms indicate 
binary zero padding (if necessary) not blanks. All characters including 
blanks, are significant in the string. 

The apostrophe may be included in an apostrophe-delimited Hollerith con- 
stant by using two consecutive apostrophes in the string. 



EXAMPLES: 



1HA 



2HA 



'ABCDE' 



Storage 



A 


b 


b 


b 


b 


b 


b 


b 


(ASCII) 


b 


A 


b 


b 


b 


b 


b 


b 


(ASCII) 


A 


B 


C 


D 


E 


b 


b 


b 



(ASCII) 
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Storage 



'ABC'DEF' 



IRA 



5LABCDE 



A 


B 


C 


t 


D 


E 


F 


b 


(ASCII) 


0000 


0000 


0000 


0041 


(HEX) 


4142 


4344 


4500 


0000 



(HEX) 



2.4,7 Hexadecimal Constants 

Hexadecimal constants are written as follows: 

X'DiD2... ! 

where the D. collectively represent a string of up to 16 hexadecimal digits 
(0-9, A-F) . Internally, hexadecimal constants are stored right justified 
in one word. They are treated as integers within expressions. 

EXAMPLES : 

X'FF 1 
X'OOOFFOOO' 
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SECTION III - SPECIFICATION STATEMENTS 

1. Introduction 

Declarations are used to supply descriptive information about the program 
other than to specify computation or other action. This descriptive 
information primarily concerns the interpretation of source program 
identifiers and object program storage allocation. 

The following declaration statements must all appear in the program prior 
to any non-declarative statements and statement function definition state- 
ments. 

Explicit type statements 
IMPLICIT statement 
DIMENSION statement 
EXTERNAL statement 
COMMON statement 
EQUIVALENCE statement 
MAXREGS statement 

DATA statements are not required to appear in the program prior to any 
non-declarative statements. 

3.1 Array Declarators 

Arrays may be defined by several statements; 

Explicit type statements 
COMMON statement 
DIMENSION statement 
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Array declarators are used in these statements. Each array declarator 
gives the array identifier and the maximum values each of its subscripts 
may assume, thus: 

identifier (maxi, max2, maxs) 

The maxima must be integers. An array may have one, two, or three dimen- 
sions. For example, the statement 

DIMENSION EDGE (10,8) 

specifies EDGE to be a two dimensional array whose first subscript may vary 
from 1 to 10 inclusive, and the second from 1 to 8 inclusive. 

EXAMPLES: 

DIMENSION PLACE (3,3,3), HI (2,4), K(256) 

Arrays may also be declared in the COMMON and Explicit type statements in 
the same way: 

COMMON X (10,4), Y, Z 
INTEGER A (7,32), B 
DOUBLE PRECISION K (6,10) 

3.1.1 Variable Dimensions 

Within a subprogram, array declarations may use integer variables provided 
that the array name and variable dimensions are formal parameters of the 
subprogram. The actual array name and values for the dummy variables are 
given by the calling program when the subprogram is called. 

Quantities needed for reference to the array are evaluated upon entry to 
the subprogram and unchanged by any subsequent modifications of the dimen- 
sion variables. 
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EXAMPLE: 

DIMENSION BETA (L,M,G) 

The identifiers BETA, L, M, and G must all be formal parameters. 

3.2 DIMENSION Statement 

The DIMENSION statement is used to declare identifiers to be array identifiers 
and to specify the number and bounds of the array subscripts. The informa- 
tion supplied in a DIMENSION statment is required for the allocation of 
memory for arrays. Any number of arrays may be declared in a single DIMEN- 
SION statement. 

FORM: 

DIMENSION Si,S 2 ,...,S k 

where S is an array declarator. 

Each array variable appearing in the program must represent an element of 
an array declared in a DIMENSION statement, unless the dimension information 
is given in another statement. When the dimension information is provided 
in a COMMON or Explicit type statement, it may not appear in a DIMENSION 
statement. The DIMENSION keyword may be abbreviated DIM. 

3.3 TYPE Declarations 

The Explicit type statements INTEGER, REAL, DOUBLE PRECISION, COMPLEX, and 
LOGICAL and the IMPLICIT statement are used to specify the type of the 
identifiers appearing in the program. 

3.3.1 Explicit Type Statements 

The general form of the Explicit type statement is: 
type ident if ieri, identifier, . . . , identifier 



n 
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Type may be INTEGER, REAL, DOUBLE, DOUBLE PRECISION, COMPLEX, or LOGICAL, 
and the identifier, are identifiers or array descriptors. 

An identifier may appear in only one Explicit type statement. Explicit 
type statements may be used to declare arrays that are not dimensioned 
in DIMENSION or COMMON statements. 



3.3.1.1 INTEGER Statement 

FORM: 

INTEGER identifier, identifier, . . . 

This statement declares the listed identifiers to be integer type with 
each datum occupying one word. 

EXAMPLES : 

INTEGER ALPHA, $PVAL 
INTEGER TABLESIZE (10) 

3.3.1.2 REAL Statement 

FORM: 

REAL identifier, identifier, .. . 

This statement declares the listed identifiers to be real type with each 
datum occupying one word in floating-point format. 

EXAMPLES : 

REAL $LOGX, MASS (10,4) 
REAL I, J, K 
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3.3.1.3 DOUBLE PRECISION Statement 

FORM: 

DOUBLE PRECISION identifier, identif ier , . . . 

This statement declares the listed identifiers to be of double precision 
type. Each datum occupies two words in floating-point format. The key- 
word PRECISION may be omitted. 

EXAMPLES : 

DOUBLE PRECISION RATE, Y, FLOW 
DOUBLE $TIME (27,9) 

3.3.1.4 COMPLEX Statement 

FORM: 

COMPLEX identifier, identifier, ... 

This statement declares the identifiers to be of complex type. Each datum 
occupies two words, two floating-point numbers representing the real and 
imaginary parts. 

EXAMPLE: 

COMPLEX ZETA, W, ROOT 

3.3.1.5 LOGICAL Statement 

FORM: 

LOGICAL identif ier, identif ier, . . . 

This statement declares the listed identifiers to be of logical type. Each 
datum occupies one word where zero represents false and non-zero represents 
true. 
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EXAMPLE: 

LOGICAL BOOL, $P, $Q, ANSWER 

3.3.2 IMPLICIT Statement 

The IMPLICIT statement defines the data identifier first-letter parameters 
for implicitly specifying data types should the programmer desire para- 
meters different from the compiler defaults, which are: 

DEFAULTS : 

Identifiers beginning with the letters I, J, K, L, M 
and N imply INTEGER data. 

Identifiers beginning with all other letters imply REAL data. 

The format of the IMPLICIT statement is: 

FORM: 

IMPLICIT type l , (Ai , A2 ,...)» type2 (Aa-Ai* ,...),... 



where: 



type, is the data type (INTEGER, REAL, DOUBLE, DOUBLE PRECISION, 
COMPLEX, LOGICAL) implied by data identifiers beginning with 
the letters listed subsequently. 

(Ai,A2,...) or (As-A^,...) are lists or ranges of alphabetic 
characters for the preceding data type. Lists of single 
characters are separated by commas; ranges of characters 
are denoted by the first and last character or the range 
(in alphabetic sequence) separated by a minus sign (e.g., 
A-D). 
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This statement causes any variable not mentioned in an Explicit type state- 
ment and whose first character is one of those listed to be typed according 
to the type appearing before the list in which the character appears. For 
an asynchronous identifier, the first character after the dollar sign is 
used. 

EXAMPLE: 

IMPLICIT INTEGER (A-C,X), DOUBLE PRECISION (D) , LOGICAL (L) 

This statement would cause the following implicit declarations to be in 
effect: 

1. Identifiers beginning with A, B, C, I, J, K, M, N, X are integer. 

2. Identifiers beginning with D are double precision. 

3. Identifiers beginning with L are logical. 

4. Identifiers beginning with E, F, G, H, 0, P, Q, R, S, T, U, V, 
W, Y, Z are real. 

3.4 COMMON Statement 

FORM: 

COMMON block-list 

The COMMON statement specifies that certain variables or arrays are to be 
stored in an area also available to other programs. By means of COMMON 
statements, a program and its subprograms may share a common storage area. 
This area is located in Data memory. 
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The common area may be divided into separate blocks identified by block 
names. A block is specified thus: 

/identifier/identifier, identifier, . . . , identifier. 

The identifier enclosed in slashes is the block name. The identifiers 
which follow are the names of the variables or arrays assigned to the 
block. These elements are placed in the block in the order- in which 
they appear in the block specification. 

The block list of the COMMON statement consists of a sequence of one or 
more block specifications. For example the statement 

COMMON/R/X,Y,T/C/U,V,W,Z 

indicates that the elements X, Y, and T, in that order, are to be placed 
in block R and that U, V, W, Z are to be placed in block C. 

Block entries concatenate throughout the program, beginning with the first 
COMMON statement. For example the statements 

COMMON/D/ALPHA/R/A,B/C/S 
COMMON/C/X,Y/R/U,V,W 

have the same effect as the statement 

COMMON/D/ALPHA/R/A,B,U,V,W/C/S,X,Y 

One block of common storage may be left unlabeled and is called blank 
common. Blank common is indicated by two consecutive slashes. For 
instance 

C0MM0N/R/X,Y//B,C,D 

indicates that B, C, and D are placed in blank common. 
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The slashes may be omitted when blank common is the first block of the 
statement. 

COMMON B,C,D 

Storage allocation for blocks with the same name begins at the same loca- 
tion for all programs executed together. For example if a program contains 

COMMON A,B/R/X,Y,Z 

as its first COMMON statement, and a subprogram has 

COMMON/R/U,V,W//D,E,F 

as its first COMMON statement, then the quantities represented by X and U 
are stored in the same location. A similar correspondence holds for A and 
D in blank common. 

Labeled blocks of a given name must have the same length in all programs 
executed together. 

Blank common may be any length in any program. 

Array names appearing in COMMON statements may have dimension information 
appended, as in a DIMENSION statement. For example 

COMMON ALPHA, T(15, 10, 5), GAMMA 

specifies the dimensions of the array T while entering T in blank common. 
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3.5 EQUIVALENCE Statement 

The EQUIVALENCE statement allows more than one identifier to represent the 
same quantity. 

FORM: 

EQUIVALENCE (Ri ,R 2 , . . . ) , (R n »R n+ i ,...),... 

where R. is a reference. (EQUIVALENCE may be abbreviated EQUIV) 

The references of an EQUIVALENCE statement may be simple variables or array 
identifiers or array element references. The subscripts of an array element 
must be integer constants. The number of subscripts must be equal to the 
array dimension or must be one. Synchronous and asynchronous variables 
may not be made equivalent to each other. 

EXAMPLE: 

EQUIVALENCE (A, B,C (3)), (T(4) ,S(1,1,2)) 

The inclusion of two or more references in a parenthesis pair indicates 
that the quantities referenced are to share same memory locations. For 
example : 

EQUIVALENCE (RED , BLUE) 
specifies that the quantities RED and BLUE are stored in the same place. 
When no array subscript is given, it is taken to be 1, thus 

EQUIVALENCE (X,Y) 

is the same as 

EQUIVALENCE (X,Y(1)) 
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Elements of multiple dimensioned arrays may be referenced with a single 
subscript by use of the element successor function. For example in the 
three-dimensional array specified by 

DIMENSION ALPHA (Ni,N 2 ,N 3 ) 

the position of element ALPHA (Ki,K2,K3> is given by element position 
« (K 3 -l)*N 1 *N 2 +(K2-l)*Ni+Ki 

Thus the sequence 

DIMENSION BETA(4),ALPHA(2,3,4) 
EQUIVALENCE (BETA(2) , ALPHA(8) ) 

specifies that BETA(2) and ALPHA (.2, 1,2) are stored in the same place. 

Since the entire arrays are shifted to satisfy the equivalence only the 
relative positions of the references are important. In the examples below 

EQUIVALENCE (BETA(l) ,ALPHA(7)) 

or 
EQUIVALENCE (BETA , ALPHA ( 7 ) ) 

will do as well. 

Note that the relation of equivalence is transitive, e.g., the two statements 

EQUI VALENCE ( A, B) , (B,C) 
EQUIVALENCE (A, B,C) 

have the same effect. 
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3.5.1 EQUIVALENCE and COMMON 

Identifiers may appear in both COMMON and EQUIVALENCE statements provided 
the following rules are observed. 

No two quantities in common may be set equivalent to one another. . 

Quantities placed in a common block by means of equivalences may cause the 
end of the common block to be extended. For example, the statements 

COMMON/R/X,Y,Z 
DIMENSION A (4) 
EQUIVALENCE(A,Y) 

causes the common block R to extend from X to A(4) , arranged as follows : 



X 

Y A(l) 

Z A(2) 

A(3) 

A(4) 



Equivalence statements which cause extension of the start of a common block 
are not allowed. For example the sequence 

COMMON/R/X,Y,Z 
DIMENSION A(4) 
EQUIVALENCE (X,A(3)) 
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is not permitted since it required block R to be arranged 

A(l) 

A(2) 
X A(3) 
Y A(4) 
X 

A(l) and A(2) extend the start of block R. 

3.6 EXTERNAL Statement 

FORM: 

EXTERNAL identifier , identifier, .... I'.^entif ier 

This statement declares the listed identifiers to be subprogram names. Any 
subprogram name given as an argument to another subprogram must appear in 
an EXTERNAL declaration in the calling program. 

EXAMPLE: 

EXTERNAL SIN, COS 



CALL TRIGF(SIN, 1.5, ANSWER) 



CALL TRIGF (COS , . 87 , ANSWER) 



END 

SUBROUTINE TRIGF (FUNC , ARG , ANSWER) 



ANSWER=FUNC(ARG) 



RETURN 
END 
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3.7 DATA Specification 

The data specification statements DATA, CONST, and BLOCK DATA are used to 
declare constants and specify initial values for variables. These values 
are compiled into the object program. They become the values assumed by 
the variables when execution begins. 

3.7.1 DATA Statement 

The data to be compiled into the object program is specified in a DATA 
statement. 

FORM: 

DATA v/d/,v/d/,... 

where v is a variable list and d is a data list. 

The variable lists in a DATA statement consist of simple variable, array 
names, or array elements separated by commas. 

Variables in common may appear on the lists only if the DATA statement 
occurs in a BLOCK DATA subprogram. Asynchronous variables may appear in 
a DATA statement. 

The data items of each data list correspond one-to-one with the variables 
of each variable list. Each data item specifies the value given to its 
corresponding variable. 

Data items may be numerical constants or alphanumeric strings. For example 

DATA ALPHA, BETA/ 5, 16. E-2/ 

specifies the value 5 for ALPHA and the value .16 for BETA. 
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Any data item may be preceded by an integer followed by an asterisk. The 
integer indicates the number of times the item is to be repeated. For 
example 

DATA A(l) ,A(2) ,A(3) ,A(4) ,A(5)/61E2,4*32E1/ 

specifies 5 values for the array A; the value 6100 for A(l) and the value 
320 for A(2) through A(5) . 

When an unsubscripted array name is included in the variable list it implies 
that all elements of the array are to be initialized. It is equivalent to 
writing out all elements of the array in sequence. 

EXAMPLE: 

The above example could be written: 
DIMENSION A(5) 
DATA A/61E2,4*32E1/ 

The form of the constant, rather than the type of the variable, determines 
the data type of the stored constant. 

Hollerith constants, or alphanumeric strings, are treated specially in the 
DATA statement. A single Hollerith constant may initialize more than one 
variable element, whereas any other type of constant corresponds to a 
single element. A Hollerith constant may not contain excess characters, 
but it will be extended with blanks so as to fill an integral number of 
elements . 

EXAMPLES: 

DIMENSION 1(2) 

DATA I/ f AVB7 

implies 1(1) = ' kWWW 

K2) - 'mmw 
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DIMENSION 1(2) 
DATA I/'ABCDEFGHIJK'/ 
implies 1(1) = 'ABCDEFGH' 
1(2) - ? IJKiMW 

DATA J/'ABCDEFGHIJKV 
is an error 

DIMENSION R(2) 
DATA R/2* f AV 
implies R(l) « ' kWUttV 
R(2) - 'kUUUV 

3.7.2 CONST Statement 

The CONST statement has the same form and general meaning as the DATA state- 
ment. 

FORM: 

CONST v/d/,v/d/,... 

where v is a variable list and d is a data li*'t. The syntax and semantics 
of these lists are identical with those of the DATA statement with one 
exception: The CONST statement causes the variables to be allocated in 
constant memory. This implies that they are not variables but symbolic 
constants. Constants may not appear in I/O statements, be passed as 
parameters, be assigned to, or be asynchronous. 

3.7.3 BLOCK DATA Statement 

FORM: 

BLOCK DATA 
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This statement declares the program which follows to be a data specification 
subprogram. Data specification for variables in common blocks requires the 
use of a BLOCK DATA subprogram. 

The first statement of the subprogram must be the BLOCK DATA statement. 
The subprogram may contain only declarative statements associated with 
the data being defined. 

EXAMPLE: 

BLOCK DATA 

COMMON/R/X,Y/C/Z,W,V 

DIMENSION Y(3) 

COMPLEX Z 

DOUBLE PRECISION X 

DATA Y(1),Y(2),Y(3),/1E-1,2*3E2/ 

DATA X,Z/11,877D0, (-1.41421,1.41421)/ 

END 

Data may be entered into more than one block of common in one subprogram, 
however, any common block mentioned must be listed in full. In the example 
above, W and V are listed in block C although no data values are defined 
for them. 
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SECTION IV - EXPRESSIONS 

4 • Introduction 

Expressions are strings of operands separated by operators, 
There are two types of FORTRAN expressions: arithmetic 
and logical. 

4 . 1 Arithmetic Expressions 

An arithmetic expression is a sequence of basic elements 
separated by arithmetic operators and parentheses in 
accordance with mathematical convention and the rules 
given below. An arithmetic expression has a single nume- 
rical value, which is the result of the calculations speci- 
fied by the quantities and operators comprising the 
expression. 

The arithmetic operators are: 



SYMBOL 


OPERATION 


+ 


Addition 


- 


Subtraction 


* 


Multiplication 


/ 


Division 


rt* 


Exponentiation 



variable, or function reference): 

2.71828 
Z(N) 

TAN(THETA) 

Compound expressions may be formed by using operators to 
combine basic elements: 
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X+3 

TOTAL/POINTS 

TAN(PI*M) 

Any expression may be enclosed in parentheses and treated 
as a basic element: 

(X+Y)/2 

(ZETA) 

COS(SIN(PI*M)+X) 

Any expression may be preceded by a + or - sign. 
For example: 

+X 

-(ALPHA BETA) 

-SQRT (-GAMMA) 

However, note that two operators must not occur consecutively 
The expression below is improper: 

X*-Y 

Use of parentheses yields the correct form: 

X*(-Y) 

By applying the rules above, all permissable arithmetic 
expressions may be formed: 

4.1.1 Operator Precedence 

If the precedence of operations is not given explicitly 
by parentheses, it is understood to be the following 
(in order of decreasing precedence) : 
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OPERATOR OPERATION 

** Exponentiation 

*and/ Multiplication and division 

+and- Addition and subtraction or negation 

For example, the expression 

A*B+C/D**E 
is executed as 

(A*B)+(C/(D**E) 

Sequences of operations of equal precedence are performed 
left to right except for exponentiation which is performed 
right to left. 

EXAMPLES 

W*X/Y/7 is evaluated as ((W*X)/Y)/7 

A**B**C is evaluated as A**(B**C) 

4.1.2 Evaluation of Mixed-Type Expr es sions 

The value of an arithmetic expression may be of integer, real, 
double precision, or complex type. The type of the expression 
is determined by the types of its elements according to the 
rules which follow. 

The arithmetic types are ranked as follows: 



RANK 


TYPE 


1 


Integer 


2 


Real 


3 


Double 


4 


Complex 



-35- 



The type of an expression is the type of the highest 
ranking element in the expression. 

Each operation within an expression is evaluated in the 
type of the highest ranking operand. Thus the evaluation 
of an expression is not changed to a higher rank until 
necessary . 

EXAMPLE 
l/J+R*DP*C 

is evaluated as 

CMPLX(FLOAT(I/J) , 0. )+CMPLX(REAL (DBLE (R) *DP) ,0.)*C 

Integer expressions are evaluated using binary integer 
arithmetic throughout. In integer arithmetic, fractional 
parts arising in division are truncated, not rounded. For 
example : 

7/3 yields 2 
3/7 yields 

All other calculations use binary floating-point arithmetic 

Conversions to higher rank are performed as follows: 

1. An integer quantity becomes the integer part of a real 
quantity. The fractional part is zero. 

2. A real quantity becomes the most significant part of a 
double precision real quantity. The least significant 
part is zero. 

3. A real quantity becomes the real part of a complex 
quantity. The imaginary part is zero. 

4. A double precision quantity is converted to single 
precision and becomes the real part of a complex 
quantity. The imaginary part is zero. 
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4.2. Logical Expressions 

There are four basic elements used FORTRAN logical expres- 
sions: logical constants, logical-type variables, logical- 
type functional references, and relations. All of these 
basic elements represent logical quantities. 

A logical quantity may have either of two values,: true 
or false. Logical quantities occupy one word of memory. 

4.2.1 Relations 

Relations are constructed from numerical expressions of 
intetger, real, or double precision type through the use 
of relational operators. The relational operators are: 

OPERATOR RELATION 

. GT . greater than 

.GE. greater than or equal to 

.LT. less than 

.LE. less than or equal to 

.EQ. equal to 

.NE. not equal to 

The enclosing periods are part of the operator and must 
be present. 

Two expressions of integer, real, or double precision type 
separated by a relational operator form a relation. For example 

X+2.LE.3*Y 
is a relation. The entire relation constitutes a basic 
logical element. 

The value of such an element is true if the relation ex- 
pressed is true a,nd false otherwise. In the example above 
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the element has the value true if X is 2 and Y is 2, and 
the value false if X is 2 and Y is 1. 

Complex operands are allowed for the opertors 
.EQ. and .NE. only. 

Relational operators have lower precedence than arithmetic 
operators. 

4.2.2. Logical Operators 

The logical operators are .NOT.,. AND. and . OR. , denoting , 
respectively, logical negation, logical multiplication 
and logical addition. The enclosing periods are part 
of the operators and must be present. The logical 
operators are defined as follows (where P and Q are 
logical expressions) ; 

.NOT.P true if P is false, 

false if P is true 
P.AND.Q true if P and Q are both true, 

otherwise false 
P.OR.Q false if P and Q are both 

false, otherwise true 

A logical expression may consist of a single logical element 
For example: 

.TRUE. 

BOOL(N) 

X.GE.3.14159 

Single elements may be combined through use of the logical 
operators .AND. and . OR. to form compound expressions, such as: 

TVAL. AND. INDEX 
BOOL(M) .OR. K.EQ. LIMIT 
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Any logical expression may be preceded by the operator. 
.NOT. as in: 

.NOT.T 

.N0T.X+7.GT.Y+Z 

BOOL(K) .AND. .NOT. (TVAL.OR.R) 

By repeated use of these rules all permissible logical 
expressions may be formed. 

When the precedence of operation is not given by parentheses, 
it is understood to be the following (in decreasing order 
of precedence) : 

OPERATOR OPERATION 

.NOT. logical negation 

.AND. logical multiplication 

.OR. logical addition 

Thus the expression 

T.AND. .NOT.S .OR. .NOT. P. AND. R 
is interpreted 

(T.AND. (.NOT.S)) .OR. ((.NOT.P) .AND.R) 

4.3. Summary of Operator Precedence 

When the precedence of operators is not given explicitly 
by parentheses, it is understood to be as follows (in order 
of decreasing precedence) : 
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OPERATOR OPERATION 

** exponential 

*,/ multiply, ' divide 

+ ,- add, subtract, negate 

• GT . , * GE • , • LT . , 

.LE. , . EQ. , .NE. relational 

.NOT. logical negation 

.AND. logical multiply 

.OR. logical add 

For example, the logical expression 

. NOT. ZETA**2+Y*MASS.GT.K-2. OR. PARITY. AND. X.EQ.Y is interpreted 
(.NOT. (((ZETA**2)+(Y*MASS)) . GT . (K-2))) .OR. (PARITY. AND. (X.EQ.Y)) 
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SECTION V - ASSIGNMENT STATEMENTS 

5 . Introduction 

Assignment statements are the basic executable statements 
of the FORTRAN language. Two types of assignment state- 
ments are available: arithmetic and Logical. 

5.1. Arithmetic Assignment Statement 

The arithmetic assignment statement specifies an arithmetic 
expression to be evaluated and a variable to which the 
expression value is to be assigned. 

FORM: 

variable = expression 

The character "=" is an operational symbol signifying 

replacement, not equality. Thus the first example below 

means "take the current value of Y, double it, and assign 
the result to Y." 

EXAMPLES: 
Y=2*Y 
A= -A 
X(N)=N*ZETA (ALPHA) *(M/PI)+(1. 0,-1.0) 

Type conversion is provided if the variable is of a type 
different from the expression. 

5 . 2 Mixed-Type Assignment 

Mixed-type assignment involves arithmetic assignment state- 
ments in which the type of the expression on the right differs 
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from the type of the variable on the left. The evaluated 
expression is converted to the type of the variable. Table 
5.1 gives the conversions for all combinations of expressions 
and variables. 



Type of Variable Type of Expression 



Conversion 



INTEGER 
INTEGER 



INTEGER 
INTEGER 

REAL 

REAL 
REAL 



COMPLEX 



COMPLEX 



(continued) 



INTEGER 
REAL 



DOUBLE PRECISION 



COMPLEX 



INTEGER 



REAL 

DOUBLE PRECISION 



REAL COMPLEX 

DOUBLE PRECISION INTEGER 

DOUBLE PRECISION COMPLEX 



INTEGER 



REAL 



Table 5.1 



None . 

Truncate fractional 
part, convert to integer 
form. 

Truncate fractional 
part, convert to integer 
form. 

Truncate fractional part 
of real part, convert to 
integer form. 
Convert to floating-point 
form. 
None. 

Take first word of ex- 
tended floating-point. 
Take real part. 
Convert to extended 
floating-point form. 
Extend real part to 
extended floating-point 
form. 

Convert to real form, 
assign to real part, 
assign zero imaginary. 
Assign to real part, 
assign zero imaginary. 
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Type o f Var ia ble 

COMPLEX 



Type of Expression 
DOUBLE PRECISION 



COMPLEX 



COMPLEX 



Conversion 

Take first word of 
extended floating-point, 
assign to real part, 
assign zero imaginary. 
None . 



Table 5.1 (continued) 

Assignment Statements of the form 
variable = Hollerith constant 
are special cases. The character string represented by 
the Hollerith constant is transferred to the variable without 
type conversion. The string is stored as indicated in sec- 
tion 2.4.6. The number of characters cannot be greater than 
the number of bytes in the variable. 

EXAMPLE: 

X='AB f 
implies the eight bytes and represented by X will appear as 



t\ 


B 


blank 


blank 


blank 


blank 




blank 


blank 



5 . 3 Logical Assignment Statement 

The logical assignment statement specifies a logical expres' 
sion to be evaluated and a logical variable to which 
the expression value is to be assigned. 



FORM: 

variable= expression 
Both the variable and expression must have logical type. 
EXAMPLES: 

LOGICAL LA,LB,LC,LD 
LA=LB . AND . LC . AND . LD 
LB^.NOT.LA 
LC=A.GT.B.OR.C.EQ.D 

-43- 



SECTION VI - CONTROL STATEMENTS 

6 . Introduction 

The normal flow of a FORTRAN program is sequentially through 
the statements in the order in which they are given to the 
compiler. By means of control statements, the programmer 
may specify the flow of the program. 

6 . 1 Statement Numbers 

FORTRAN statements may be given numbers to be referenced 
by control statements. A statement number is written as an 
unsigned integer of five digits or less. Leading zeros and 
embedded blanks are ignored. 

Although statement numbers are written as integers, they do 
not represent numerical quantities. Statement numbers repre- 
sent statement labels, a distinct basic quantity. Statement 
numbers are used for program control, not numerical calculation 

Statement numbers must be unique, i.e., no two statements may 
have the same number. 

6.2 GO TO Statements 

GO TO statements unconditionally transfer control from one 
part of the program to another. There are three forms of the 
GO TO statement: unconditional, computed, and assigned. 

6.2.1 Unconditional GO TO Statement 

FORM: 

GO TO n 
where n is a statement number. 
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This statement transfers control to the statement num- 
bered n. 

EXAMPLE: 

GO TO 345 

6.2.2 Computed GO TO Statement 

FORM: 

GO TO (n x , n 2 ,... ,n K ) , variable 
where n r ,n^,,..,njj are statement number. 

The variable must be of integer type. 

The statement transfers control to the statement numbered 
n^,n£,...4£ if the variable has the value 1, 2, . . . k, re- 
spectively. Values outside the range 1 through k cause 
transfer of control to the following statement. The comma 
preceding the variable may be omitted. 

EXAMPLES: 

GO TO (22,3,7) , SWITCH 
GO TO (1,2,62,78)Y 

6.2.3 Assigned GO TO Statements 

FORMS: 

GO TO variable 

GO TO variable, (n x ,n 2 , ... ,n^) 
where rij ,n 2 , . . . n K are statement numbers. 

The second form is allowed for compatibility only; the labels 
nj,...n K are not used. The comma following the variable may 
be omitted. 

The variable must be a scalar of integer type and must not be 
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an asynchronous variable. This stsUient transfers control 
to the statement whose number was last assigned to the 
variable. The assignment must take place in a previously 
executed ASSIGN statement. 

The variable is a control variable, having a label as a 
value, not a numerical quantity. 

EXAMPLES: 

GO TO ERROR 

GO TO X(100.200) 



6.2.4 ASSIGN Statement 



FORM: 

ASSIGN statement number TO variable 

The variable must be a scalar of integer type and must 
be an asynchronous variable. This statement assigns the 
value of the variable for a subsequent assigned GO TO 
statement. The statement number represents the statement 
to which the assigned GO TO will transfer control. 

EXAMPLES: 

ASSIGN 7 TO LABEL 
ASSIGN 13 TO ERROR 

6 . 3 Arithmetic IF Statement 

FORM: 

IF (expression) nj ,n 2 ,n 3 
where n 1 ,n 2 >n 3 are statement numbers. 

This statement transfers control to the statement numbered 
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ni 112 or n3 if the value of the expression is less than, 
equal to, or greater than zero, respectively. The expres- 
sion must be of integer, real, or double precision type. 

EXAMPLES: 

IF (ETA)4,7,12 

IF (KAPPA-L(10))20,14,14 

6.4. Logical Statement 

FORM: 

IF (expression) S 
where S is a complete statement. 

The expression must be a logical expression. S may be any 
imperative (executable) statement other than a DO statement 
or another logical IF statement. 

If the value of the expression is false, then control passes 
to the next sequential statement. 

If the value of the expression is true, statement S is executed 
After execution of S, control passes to the next sequential 
statement unless S is an arithmetic i* v statement or GO TO type 
statement in which case control is transferred as indicated. 

As an example, consider the statements 

IF(B)Y*=X*SIN(Z) 
W=Y**2 

If the value of B is true the statements Y=»X*SIN(Z) and 
W=Y**2 are executed in that order. 

If the value of B is false, the statement Y=X*SIN(Z) 
is not executed. 
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EXAMPLES: 

IF(T.0R.S).X=Y+1 
IF(Z.GT.X(K))CALL SWITCH(X,Y) 
IF(K.EQ.INDEX)GO TO 15 

6.5 DO Statement 

FORMS: 

DO n index ^initial, limit 
DO n index =initial, limit, step 
where n is a statement number. 

The index must be a simple integer variable, and must not 
be an asynchronous variable. The initial and limit may 
be integer simple variables, signed integer constants, or 
integer expressions. The step must be an integer simple 
variable, an integer constant or an integer expression. 
If the step is not given, it is understood to be one. 

The DO statement causes the statements which follow, up 
to and including the statement numbered n, to be executed 
repeatedly. This group of statements is called the range 
of the DO statement. Initially, the initial value is 
assigned to the index. Thereafter, after each execution 
of the range, the step value is added to the index value 
and the result assigned to the index. 

Prior to each execution of the range, the index value is 
compared to the limit value. If the index value does 
not exceed the limit value, the range is executed. This 
differs from Standard FORTRAN in which the range is always 
executed once before the first test. 

After the last execution of the range, control passes 
to the statement immediately following the range. The 
exit from the range is called the normal exit. 
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Exit may also be accomplished by a transfer from within 
the range. 

The range of a DO statement may include other DO state- 
ment provided that the range of each contained DO state- 
maent is entirely within the range of the containing 
DO statement. 

Within the range of a DO statement, the index is available 
for use as an ordinary variable. After a transfer exit 
from the range, the index retains its current value 
and is available for use as an ordinary variable. After 
a normal exit from the range, the index retains the value 
which caused the exit. 

The values of the index, limit, and step may be altered 
within the range of the DO statement. Altering the value 
of limit or step does not affect the loop. 

The range of a DO statement must not end with a GO TO 
type statement or an arithmetic IF statement. A logical 
IF statement is allowed as the last statement of the range, 
provided the logical IF does not contain a GO TO type 
statement or an arithmetic IF. In this case, control is 
transferred thus: 

The range is considered ended when and if control 
would normally pass to the statement following the 
logical IF statement. 

As an example, consider the sequence: 

DO 5K = 1,4 

5 IF(X(K) .GT.Y(K))Y(K) = (K) 

6 . * . 
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Statement 5 is executed four times, whether the state- 
ment Y(K) = X(K) is executed or not. 

Statement 6 is not executed until statement 5 has been 
executed four times. 

EXAMPLES: 

DO 22 L » 1,30 

DO 45K = 2, LIMIT, 3 

6.6. CONTINUE Statement 

This statement is a dummy statement, used primarily as 
a target point for transfers, particualrly as the last 
statement in the range of a DO statement. For example 
in the sequence 

FORM: 

CONTINUE 

DO 7K = START, END 



IF(X(K))22,13,7 
7 CONTINUE 

a positive value X(K) begins another execution of the range 
The CONTINUE provides a target address for the IF statement 
and ends the range of the 

6. 7 PAUSE Statement 

FORMS: 

PAUSE 
PAUSE n 

PAUSE character string 
where n is an integer constant. 
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This statement causes the operand to be displayed on the system 
output device. If no operand is given, zero is displayed. 

EXAMPLE: 

PAUSE 167 

6.8 STOP Statement 

FORMS : 

STOP 
STOP n 

STOP 'character string* 
where n is an integer constant. 

This statement terminates the program. It does not stop the 
system and may be used to transfer control to the operating 
system. 

6. 9 PURGE Statement 

FORM: 

PURGE Ni N 2 , N n 

where the % are asynchronous variable, 

The statement causes the access state of all the named 
variables to be set empty; regardless of their existing state. 
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SECTION VII - INPUT/OUTPUT 

7. Introduction 

Input/output statements specify the transfer of information between com- 
puter memory and input/output devices, or between one part of computer 
memory and another. These statements also allow the program to mani- 
pulate I/O devices. 

Information may be transferred in two different forms: formatted and 
unformatted. The unformatted form involves no data conversion, data 
is transferred in its internal format. Formatted data is converted 
from internal to external form, or vice versa, under control of a 
FORMAT specification in the program. 

7.1 Input/Output Lists 

Input /output statements may contain a list of variables which are to 
receive values on input or are to provide values for output. The list 
of a transmission statement specifies the order of transmission of the 
variable values. During input, the new values of listed variables may 
be used in subscript or control expressions for variables appearing later 
in the list. For example 

READ (1,3)L,A(L),B(L+1) 

reads a new value of L and uses this value in the subscripts of A and B. 

The transmission of array variables may be controlled by indexing similar 
to the indexing in the DO statement. The list of controlled variables, 
followed by the index control, is enclosed in parentheses and the whole 
acts as a single element of the list. For example 

READ (7,23)(X(K),K=1,4) 

-52- 



is equivalent to 

READ (7,23)X(1),X(2),X(3),X(4) 
The initial, limit, and step values are given as in the DO statement: 

READ (4,2)N,(GAIN(K) ,K=1,H,N) 
The indexing may be compounded as in the following: 

READ (1,13)((MASS(K,L)K=1,5),L«1,4) 
This statement reads the elements of array MASS in the order 

MASS(1,1) ,MASS(2,1) , . . . ,MASS(5.1) ,MASS(1.2) , . . . ,MASS(5,4) 

If an entire array is to be transmitted, the indexing may be omitted and 
only the array identifier written. The array is transmitted in order of 
increasing subscripts with the first subscript varying most rapidly. 
Thus the example above can be written 

READ(1,13)MASS 

7.2 Input/Output Statement Parameters 

Many input/output statements have similar formats. The following defini- 
tions apply to all input/output statements in which they may appear: 

u logical I/O unit number which may be an unsigned integer 
constant or an integer simple variable and may not be an 
asynchronous variable. The correspondence between unit 
number and actual I/O device is determined by the system 
configuration and the operating system. 
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f format declaration identifier which is the statement number 
of a FORMAT statement in the program. 

list I/O list as defined in section 7.1. 

Si statement number to which program control is transferred in 
the event an end~of-file indication is detected on the I/O 
unit while processing the statement. - 

S2 statement number to which program control is transferred if 
any error is detected while processing the statement. 

7.3 Sequential Input/Output Statements 

These statements treat input and output to and from I/O units as if the 
units contained sequential files, each composed of an ordered set of 
records. Each time a READ or WRITE statement is executed, at least one 
record is processed. As each record is processed, the file is positioned 
to read or write the next sequential record. 

7.3.1 READ Statement 

FORMS : 

READ (u,f,END=Si,ERR=S 2 ) list 
READ (u,END=S 1 ,ERR«S 2 ) list 

where u is an I/O unit designation and f is a format reference. 

The parameters END=Si and ERR 3 ^ are optional and their order may be reversed. 

The READ statement causes information to be rsad from the I/O unit designated 
and stored in memory as values of the variables in the list. 
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In transmitting formatted data, the conversion from external to internal 
form is specified by the format referenced (first form) . 

When binary data is transmitted the format reference is omitted (second 
form) . 

EXAMPLES : 

READ (1,15 END=100)ETA,PI 

READ (K+L,10)GSIN,ZAI 

READ (M,FMT,ERR=*999,END=100) (TABL(K) ,K=1,M) 

READ (TAPE) (TEMP (L),L=1, 100) 

7.3.2 WRITE Statement 

FORMS: 

WRITE (u,f,END=Si,ERR=*S 2 ) list 
WRITE (u,END=S 1 ,ERR«S 2 ) list 

where u is an I/O unit designation and f is a format reference. 

The END=Si and ERR=S2 parameters are optional and their order may be 
reversed. 

The WRITE statement causes the values of the variables in the list to be 
transmitted from memory to the designated I/O unit. 

In transmitting formatted data, the conversion from internal to external 
form is specified by the format reference (first form) . 
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When binary data is transmitted to the format reference is omitted (second 
form) . 

EXAMPLES : 

WRITE (2,15)ZILCH 

WRITE (K3,4,)(A(K),K=2,20),ICHI 
WRITE (4)NUMB,(SYMB(J),J=1,NUMB) 

7.4 FORMAT specifications 

The format designator f appearing in formatted I/O statements is the 
statement label of a FORMAT statement. 



FORM: 



F0RMAT(Si,S 2 ,...,S k ) 



where S is a data field specification. 
7.4.1 N umeric Fields 
Conversion of numerical data may be one of six types. 

1. type - E 

internal form - binary floating-point 
external form - decimal floating-point 

2. type - F 

internal form - binary floating-point 
external form - decimal fixed-point 

3. type - G 

internal form - binary floating-point 

external form - decimal fixed-point or floating-point 
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4. type - I 

internal form - binary integer 
external form - decimal integer 

5. type - Z 

internal form - binary integer 
external form - hexadecimal integer 

These types of conversion are specified by the forms 

1. Ew.d 

2. Fw.d 

3. Gw.d 

4. Iw 

5. Zw 

respectively. The letter E, F, G, I, or Z designates the conversion type; 
w is an integer specifying the field width; d is in integer specifying the 
number of decimal places to the right of the decimal point. As an aid in 
conversion, the letter is allowed and treated as a Z. For example, the 
statement 

FORMAT (I5,F10.2,E25.15) 
could be used to output the line 

32 -17.60 5.9625478777541E03 

on the output listing. 

The type of conversion used should correspond to the type of the variable 
in the input/output list. I conversion is used for integer type variables 
and E, F, or G conversion is used for real type variables. There is no 
format conversion for double precision variables. 
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The decimal fixed-point number (type F) has a decimal point but no exponent, 
whereas the decimal floating-point (type E) has an exponent. On output, 
the exponent always has the form shown, i.e., an "E" followed by a signed, 
two digit integer. On input, however the "E" or the "+" sign, or the 
entire exponent may be omitted on the external form. For example, the 
following are all valid E15.6 fields: 

.317250+2 
.317250E2 
.042739-45 
31064 

The field width w includes all of the characters (decimal point, signs, 
blanks, etc.) which comprise the number. If a number is too long for 
its specified field, the excess characters are lost. Since numbers are 
right justified in their fields, the loss is from the most significant 
part of the number. 

During input, the appearance of a decimal point (.) in an E, or F type 
number overrrides the d specification of the field. In the absence of 
an explicit decimal point, the point is positioned d places from the 
right of the field, not counting the exponent, if present. For example, 
a number with external appearance 271828E-1 and specification E12.5 
is interpreted as 2.71828E-1. 

7.4.2 Scale Factors 

Scale factors may be specified for E, and F type conversions. A scale 
factor is written nP where P is the identifying character and n is a 
signed or unsigned integer specifying the scale factor. 

For F type conversion the scale factor specifies a power of ten such that 

external number = (internal number)* (power of ten) 
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For E type conversions, the scale factor multiplies the number by a power 
of ten but the exponent is changed accordingly, leaving the number un- 
changed except in form. For example if the statement 

FORMAT (F8.3,E16.5) 
corresponds to the line 

26.451 -4.1321E-02 
the the statement 

FORMAT (-1PF8.3,1PE16.5) 
corresponds to the line 

2.645 -41.32100E-03 

The default scale factor is 0. However, once a scale factor is given, it 
holds for all following E.and F type conversions within the same format. 
The scale factor is reset to zero by giving a scale factor of zero. 

Scale factors have no effect on I conversions. 

7.4.3 G-Fields 

Output under control of a G-field is dependent on the magnitude of the 
floating-point number being converted. Where m represents the magnitude 
of the number, the following table shows the relationship between m and 
the conversion field to be used. 
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Magnitude Conversion Field 

0.1<m<l F(w-4).d,4X 

Km<10 F(w-4) . (d-1) ,4X 



10 d ~ 2 <to<10 d ~ l F(w-4).1,4X 

10 d " 1 <m<10 d F(w-4) .0,4X 

nK.l or m>10 sEw.d 

s is the current scale factor and applies only when the E conversion field 
is used, 4X denotes a field of four spaces. 

Input under control of a G-field is the same as for the F-field. 

7.4.4 Logical Fields 

Logical data can be transmitted in a manner similar to numeric data by use 
of the form: 



Lw 



where L is the control character and w is an integer specifying the field 
width. 

Data is transmitted as the value of a logical variable in the input/output 
list. 

On input, the data field is inspected for a T or F. If one is found the 
value of the logical variable is stored as true or false, respectively. 
If the data field contains no T or F, a value of false is stored. 

On output, w-1 blanks followed by T or F is output if the value of the 
logical variable is true or false, respectively. 
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7.4.5 Alphanumeric Fields 

Alphanumeric data can be transmitted in a manner similar to numeric data 
by use of the form Aw or Rw; A and R are the control characters and w 
is the number of characters in the field. The alphanumeric characters 
are transmitted as the value of a variable in an input/output list. The 
variable may be of any type. For example, the sequence 

READ (2,5)V 
5 FORMAT (A4) 

causes four characters to be read and placed in memory as the value of the 
variable V. 

The character information is transferred as 8-bit ASCII characters, stored 
8 characters per 64-bit word. 

Although w may have any value, the number of characters transmitted is 
limited by the maximum number of characters which can be stored in the 
space allotted for the variable. This maximum depends on the variable 
type. If w exceeds the maximum, leading characters are lost on input 
and replaced with blanks on output. When w is less than the maximum, 
the A format causes left justification with blanks filled on input. 
Only the left-most w characters are used for output. The R format 
causes right justification with binary zeros filled on input. Only 
the right-most w characters are used for output. 

7.4.6 Alphanumeric Constant Fields 

An alphanumeric constant may be specified within a format by preceding 
the alphanumeric string by the form nH. H is the control character and 
n is the number of characters in the string, counting blanks. For example, 
the statement 

FORMAT (17H PROGRAM COMPLETE) 
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can be used to output 

PROGRAM COMPLETE 

on the output listing. 

Alphanumeric strings delimited by single quotes may be used in the same 
manner . 

7.4.7 Mixed Fields 

An alphanumeric format field may be placed among other fields of the format, 
For example, the statement 

FORMAT (15, 8H F0RCE=F10.5) 

can be used to output the line 

22 F0RCE=17. 68901 

Note that the separating comma may be omitted after an alphanumeric format 
field. 

7.4.8 Blank or Skip Fields 

Blanks may be introduced into an output record or characters skipped on an 
input record by use of the specification nX. The control character is X 
and n is the number of blanks or characters skipped, n must be greater 
than zero. For example, the statement 

FORMAT (5H STEPI5,10X,3HY=F7.3) 
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may be used to output the line 
STEP 28flMM*WM*Y«-3.872 
where ten blanks separate the two quantities. 

7.4.9 Tabulation 

The position in the record where the transfer of data is to begin can be 
specified by T format conversion. The specification is Tn where n is the 
character position. For printed output, the first character is for carriage 
control and should not be counted. 

EXAMPLE: 

FORMAT (T20, 'NAME 1 , T40, f AGE', T1,6H GRADE) 

would print a line: 

Position 1 Position 19 Position 39 

4' 4" 4 

GRADE NAME AGE 

7.4.10 Repetition of Field Specifications 

Repetition of a field specification may be specified by preceding the control 
character E, F, G, I by an unsigned integer giving the number of repetitions 
desired. For example 

FORMAT (2E12.4,3I5) 

is equivalent to 

FORMAT (E12.4,E12.4,I5,I5,I5) 
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7.4.11 Repetition of Groups 

A group of field specifications may be repeated by enclosing the group in 
parentheses and preceding the whole with the repetition number. For example 

FORMAT (2I8,2(EI5.5,2(F8.3))) 

is equivalent to 

FORMAT (2I8,E15.5,2F8.3,E15.5,2F8.3) 

7.4.12 Complex Fields 

Complex quantities are transmitted as two independent real quantities. The 
format specification is given as two successive real specifications or one 
repeated real specification. For instance, tY»3 statement 

FORMAT (2E15.4,2(F8.3,F8.5)) 

could be used in the transmission of three complex quantities. 

7.4.13 Multiple-Record Formats 

To handle a group of input/output records where different records have 
different field specifications, a slash (/) is used to indicate a new 
record. For example, the statement 

FORMAT (318/15, 2F8. 4) 

is equivalent to 

FORMAT (318) 
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for the first record and 

FORMAT (I5,2F8.4) 

for the second record. 

The separating comma may be omitted when a slash is used. 

Blank records may be written on output or records skipped on input by using 
consecutive slashes. 

Both the slash and the closing parenthesis at the end of the format indicate 
the termination of a record. If the list of an input/output statement 
dictates that transmission of data is to continue after the closing paren- 
thesis of the format is reached, the format is repeated from the last open 
parenthesis level of one or zero. If this parenthesis is preceded by a 
repeat specification, the repeat specification is reused. Thus the statement 

FORMAT (F7.2,2(E15.5,E15.4),I7) 
causes the format 

F7.2,2(E15.5,E15.4),I7 
to be used on the first record and the format 

2(15.5,E15.4),I7 
on succeeding records. 
As a further example, consider the following statement. 

FORMAT (F7.2/(2(E15.5,E15.4),I7)) 
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The first record has the format below. 

F7.2 
Successive records have the following format. 

2(E15.5,E15.4),I7 

7 . 4 . 14 Carriage Control for Printing 

Every record that is transmitted to a listing device for printing is assumed 
to have a carriage control character as the first character of the record. 
The carriage control character itself is not printed. The carriage control 
characters are: 

Character Function Before Printing 

blank Space one line 

Space two lines 

1 Skip to first line of next page 

any other character is treated as a blank. 

EXAMPLE: 

10 F0RMAT(9H1 PAGE ,I3/1H0) 

7.5 Auxiliary I/O Statements 

These statements are used to control the positioning and file marking of 
sequential files. 
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7.5.1 REWIND Statement 

FORM: 

REWIND u 

where u is an I/O unit designation. 

This statement directs the I/O unit designated to reposition to the first 
record, u must not be an asynchronous variable. 

EXAMPLES: 

REWIND 2 
REWIND K 

7.5.2 BACKSPACE Statement 

FORM: 

BACKSPACE u 

where u is an I/O unit designation. 

This statement directs the I/O unit designated to backspace one record, u 
must not be an ansynchronous variable. 

EXAMPLES: 

BACKSPACE 5 
BACKSPACE N 
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7.5.3 END FILE Statement 

FORM: 

END FILE u 

where u is an I/O designation. 

The statement directs the I/O unit designated to write an end-file mark, 
u must not be an asynchronous variable. 

EXAMPLE: 

END FILE 4 
END FILE T 
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SECTION VIII - PROGRAM UNITS 

8. Introduction 

A FORTRAN program consists of one main program and, optionally, SUBROUTINE 
subprograms, FUNCTION subprograms, and BLOCK DATA subprograms. Each of 
these is termed a "program unit", 

8.1 PROGRAM Statement 

FORM: 

PROGRAM identifier 

The PROGRAM statement defines the program name that is used as the entry- 
point name for the object module. The identifier must not appear any- 
where else in the program unit. This statement, if present, must be the 
first statement of a main program. If not present, the main program name 
defaults fo F%MAIN. 

8.2 END Statement 

FORM: 

END 

The END statement must be the last physical statement of each program unit 
It informs the compiler of the end of the program unit. The END statement 
must be on a single source line; continuation lines are not allowed. 

8.3 RETURN Statement 

FORM: 

RETURN 
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This statement returns control from a FUNCTION or SUBROUTINE subprogram 
to the calling program unit. Normally, the last statement executed in a 
subprogram is a RETURN statement. It need not be the last statement of 
the program. Any number of RETURN statements may be used. 

8.4 RESUME Statement 

FORM: 

RESUME 

This statement is used to place a function or subroutine in the asynchronous 
or parallel mode. It allows the calling program unit to continue execution; 
however, the function or subroutine continues to execute as well until a 
RETURN statement is executed. 

8.5 Subprogram Communications 

The main program and subprograms communicate with each other by means of 
COMMON variables and parameters. If the means of communication is by 
parameters, the arguments of the subroutine or function call are known as 
actual parameters. Corresponding arguments in the subroutine or function 
argument list are known as formal parameters. 

8.5.1 Actual Parameters 

The actual parameters which appear in a subroutine call or a function 
reference may be any of the following: 

An arithmetic expression 

A logical expression 

A constant 

A simple variable 

An array element reference 

An array name 

A FUNCTION name 

A SUBROUTINE name 
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8.5.2 Formal Parameters 

The formal parameters appearing in the parenthetical list of a FUNCTION 
or SUBROUTINE statement may be any of the following: 

An array name 

A simple variable 

A subprogram name 

(either function or subroutine) 

The formal parameters are replaced at each execution of the subprogram by 
the actual parameters supplied in the CALL statement or function reference, 

Formal parameters representing array names must appear within the sub- 
program in type or DIMENSION statements giving dimension information. In 
a type or DIMENSION statement, formal parameters may be used to specify 
variable dimensions for array name formal parameters. Variable dimensions 
may be given only for arrays which are formal parameters. 

Within a subprogram, the use of formal parameters is restricted as follows; 

1. Formal parameters must not appear in COMMON statements. 

2. Formal parameters must not appear in EQUIVALENCE statements. 

3. Formal parameters must not appear in DATA statements. 

8.5.3 Correspondence Between Actual and Formal Parameters 

When a subprogram is called, the formal parameters must agree with the 
actual parameters as to number, order, type, and length. For example, if 
an actual parameter is an integer constant, then the corresponding formal 
parameter must be of INTEGER type. 

Also, the formal and actual parameters must be either both synchronous or 
both asynchronous, they must not be mixed. 
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If a formal parameter is an array name, the corresponding actual parameter 
may be either an array name or an array element. 

If a formal parameter is assigned a value in the subprogram, the corres- 
ponding actual parameter must be a simple variable array element, or 
array name. A constant or expression should not be used as an actual 
parameter if the corresponding formal parameter may be assigned a value. 

8.6 Statement Function Definition Statement 

FORM: 

identifier (identifier, identifier, . . . )=expression 

This statement defines an internal subprogram. The entire definition is 
contained in a single statement. The first identifier is the name of the 
subprogram being defined. 

Statement function subprograms are functions; they are single-value and 
must have at least one argument. The type of the function is determined 
by the type of the function identifier. 

The identifiers enclosed in parentheses represent the arguments of the 
function. These are formal parameters which have meaning and must be 
unique only within the statement. They may be identical to identifiers 
of the same type appearing elsewhere in the program. These identifiers 
must agree in order, number, type, and length with the actual parameters 
given at execution time. 

The use of a parameter in the defining expression is specified by the use 
of its parameter identifier. Expressions are the only permissible arguments 
of internal functions; hence the parameter identifiers may appear only as 
simple variables in the defining expression. They may not appear as array 
identifier. 
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Identifiers appearing in the defining expression which do not represent 
parameters are treated as ordinary variables. 

The defining expression may include references to external functions or 
other previously defined internal functions. 

All statement function definition statements must precede the first exe- 
cutable statement of the program. 

EXAMPLES : 

SSQR(K)=K*(K+l)*(2K+l)/6 

NOR(T,S)«.NOT. (T.OR.S) 

ACOSH (X) » (EXP (X/ A)+EXP (-X/A) ) / 2 

In the last example above, X is a parameter identifier and A is an ordinary 
identifier. At execution, the function is evaluated using the current 
value of the quantity represented by A. 

8.7 FUNCTION Subprogram s 

A FUNCTION subprogram is a function; it returns a single value and is 
referenced as a basic element in an expression. A FUNCTION subprogram 
begins with a FUNCTION declaration and returns control to the calling 
program by means of a RETURN or RESUME statement. It is a program unit 
and, consequently, must terminate with an END statement. 

8.7.1 FUNCTION Statement 

FORM: 

FUNCTION identifier (identifier, identifier, . . .) 

This statement declares the program which follows to be a function subpro- 
gram. The first identifier is the name of the function being defined. This 
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identifier must appear as a simple variable srvl be assigned a value during 
execution of the subprogram. This value is the function value. 

Identifiers appearing in the list enclosed in parentheses are formal 
parmeters representing the function arguments. 

EXAMPLE: 

FUNCTION FLOAT (I) 

FLOAT=I 

RETURN 

END 

8.7.2 FUNCTION Type 

The type of the function is the type of identifier used to name the function. 
This identifier may be typed implicitly or explicitly in the same way as 
any other identifier. Alternately, the function may be explicitly typed 
in the FUNCTION statement itself by replacing the word FUNCTION with one 
of the following: 

INTEGER FUNCTION 

REAL FUNCTION 

DOUBLE FUNCTION 

DOUBLE PRECISION FUNCTION 

COMPLEX FUNCTION 

LOGICAL FUNCTION 

for example, the statement 

COMPLEX FUNCTION HPRIME(S,N) 

is equivalent to the statements 

FUNCTION HPRIME(S,N) 
COMPLEX HPRIME 
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EXAMPLES : 

FUNCTION MAY (RANG, XP,YP,ZP) 
REAL FUNCTION COT(ARG) 

8.7.3 Library Functions 

The FORTRAN system supplies a library of standard functions which may be 
referenced from any program. Appendix A lists these library functions. 
These are divided into two sets: basic external functions and intrinsic 
functions. The basic external functions are called by the object program 
in the same manner as normal, user-supplied functions. Intrinsic func- 
tion names are known to the compiler and intrinsic function references 
may be treated in non-standard ways (such as expanding the function 
in-line) . The programmer can supply his own function in place of an 
intrinsic function by including the name in EXTERNAL statements in all 
calling programs. 

8.8 SUBROUTINE Subprograms 

A SUBROUTINE subprogram is not a function; it can be referred to only 
by a CALL or CREATE statement. A SUBROUTINE subprogram begins with a 
SUBROUTINE declaration and returns control to the calling program by 
means of a RETURN or RESUME statement. 

8.8.1 SUBROUTINE Statement 

FORMS: 

SUBROUTINE identifier 

SUBROUTINE identifier (identifier , identifier , . . . ) 

This statement declares the program which follows to be a SUBROUTINE sub- 
program. The first identifier is the subroutine name. The identifiers 
in the list enclosed in parentheses are formal parameters. 
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A SUBROUTINE subprogram may use one or more of its formal parameters to 
represent results. The subprogram name is not used for return of results. 

A SUBROUTINE subprogram need not have any parameters at all. 

EXAMPLES : 

SUBROUTINE EXIT 

SUBROUTINE FACTOR (CEF,N, ROOTS) 

SUBROUTINE RESIDUE (NUM,D,DEN,M,RES) 

8.8.2 CALL Statement 

FORMS: 

CALL identifier 

CALL identifier (argument, argument, . . . , argument) 

The CALL statement is used to transfer control to a subroutine subprogram. 
The identifier is the subprogram name. 

The parameters may be expressions, array identifiers, alphanumeric strings, 
or subprogram identifiers, as in the case of a function reference. Unlike 
a function, however, a subroutine cannot be referenced as a basic element 
in an expression. A subroutine may use one or more of its arguments to 
return results to the calling program. If no arguments at all are required, 
the first form is used. 

EXAMPLES : 

CALL EXIT 

CALL SWITCH (SIN,2.LE.BETA,X**4,Y) 

CALL MULT (A,B,C) 

The identifier used to name the subroutine is not assigned a type and has 
no relation to the types of the arguments. 
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8.8.3 CREATE Statement 



FORMS: 



CREATE identifier 

CREATE identifier (argument, argument, . . . , argument) 

The CREATE statement is used to execute a subroutine as a parallel pro- 
cess. The identifier is the subroutine name. 

The CREATE statement is the only way to directly initiate execution of 
an asynchronous subroutine. Note, however, that parallel execution is 
also obtained through the use of the RESUME statement. 

EXAMPLES : 

CREATE GRAPH ($IN,X) 
CREATE PROC 
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APPENDIX A 
LIBRARY FUNCTIONS 



A.l Intrinsic Functions 



FUNCTION 




Type of 
Parameter 


ABS(a) 




Real 


IABS(a) 




Integer 


DABS (a) 




Double 


AINT(a) 




Real 


INT(a) 




Real 


IDINT(a) 




Double 


AM0D(ai,a 2 ) 




Real 


M0D(ai,a 2 ) 




Integer 


AMAX0(ai,a 2 ,. 


..) 


Integer 


AMAXl(ai,a 2 ,. 


..) 


Real 


MAX0(ai,a 2 ,.. 


.) 


Integer 


MAXl(ai ,a 2 , . . 


.) 


Real 


DMAXl(a!,a 2 ,. 


..) 


Double 


AMIN0(ai,a 2 ,. 


..) 


Integer 


AMINl(ai,a 2 , . 


..) 


Real 


MINO(ai,a 2 ,.. 


.) 


Integer 


MINl(ai,a 2 ,.. 


.) 


Real 


DMINl(ai,a 2 ,. 


..) 


Double 



Type of 
Result 

Real 

Integer 

Double 

Real 

Integer 

Integer 

Real 
Integer 

Real 

Real 

Integer 

Integer 

Double 

Read 

Real 

Integer 

Integer 

Double 



Definition 



Truncation 



ai(mod a 2 ) 



Max(ai ,a 2 , . . .) 



Min(ai,a 2 , . . .) 



FLOAT (a) 



Integer 



Real 



Conversion from integer 
to real 



FULL(a) 



Any 

asynchronous 

type 



Logical Test FULL access state 
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Function 



EMPTY (a) 



Type of 
Parameter 

Any 

asynchronous 

type 



Type of 

Result Definition 

Logical Test EMPTY access state 



IFIX(a) 



Real 



Integer Conversion from real to 
integer 



SIGN(ai,a 2 ) 


Real 


Real 


Sign of a 2 times |ai 


ISIGN(ai,a2) 


Integer 


Integer 




DSIGN(a!,a 2 ) 


Double 


Double 




DIM(ai,a 2 ) 


Real 


Real 


ai-Min(ai,a 2 ) 


IDIM(ai,a 2 ) 


Integer 


Integer 





SNGL(a) 



Double 



Real 



Conversion from double 
to real 



REAL (a) 
AIMAG(a) 



Complex 
Complex 



Real 



Real 



Obtain real part of complex 



Obtain imaginary part of 
complex 



DBLE(a) 

CMPLX(ai,a 2 ) 
CONJG(a) 



Real 



Real 



Complex 



Double Conversion from real to 
double 



Complex ai + a 2 /^X 



Complex Obtain conjugate of 
complex 



I0R(ai,a 2 ) 
LAND(ai,a 2 ) 



Integer 
Integer 



Integer Inclusive OR 



Integer Logical AND 
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Function 
NOT(ai) 

IEOR(ai,a 2 ) 

ISHFT(ai,a 2 ) 



Type of 
Parameter 

Integer 
Integer 
Integer 



Type of 
Result 



Definition 



Integer Logical negation 



Integer 



Exclusive OR 



Integer Shift ai by a 2 bits 
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A. 2 Basic External Functions 



Function 


Type of 
Parameter 


Type of 
Result 


Definition 


EXP (a) 


Real 


Real 


a 
e 


DEXP(a) 


Double 


Double 




CEXP(a) 


Complex 


Complex 




ALOG(a) 


Real 


Real 


In (a) 


DLOG(a) 


Double 


Double 




CLOG(a) 


Complex 


Complex 




AL0G10(a) 


Real 


Real 


logio(a) 


DLOGlO(a) 


Double 


Double 




SIN(a) 


Real 


Real 


sin(a) 


DSIN(a) 


Double 


Double 




CSIN(a) 


Complex 


Complex 




COS (a) 


Real 


Real 


cos (a) 


DCOS(a) 


Double 


Double 




CCOS(a) 


Complex 


Complex 




TANH(a) 


Real 


Real 


tanh(a) 


SQRT(a) 


Real 


Real 


(a) l/2 


DSQRT(a) 


Double 


Double 




CSQRT(a) 


Complex 


Complex 




ATAN(a) 


Real 


Real 


arctan(a) 


DATAN(a) 


Double 


Double 




ATAN2(ai,a 2 ) 


Real 


Real 


arctan(ai/a2) 


DATAN2(ai,a 2 ) 


Double 


Double 
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Function 
DM0D(ai,a 2 ) 

CABS (a) 



Type of 
Parameter 

Double 



Complex 



Type of 
Result 

Double 



Definition 



ai(mod &z) 



Complex Absolute value of Complex 
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APPENDIX B 
FORTRAN LISTING FORMAT 

B.l Introduction 

The printed output from the FORTRAN compiler is organized as follows for 
each program or subprogram in a compilation: 

1. Program listing. 

2. Allocation of variables for each COMMON block and for non-COMMON 
variables - gives location within the block, variable name, data 
type, size, and whether scalar or array. 

3. Cross-reference (optional)- lists each symbol followed by the line 
number of each statement in which the symbol appears. A line 
number surrounded by slashes indicates appearance in a speci- 
fication statement. A line number surrounded by asterisks 
indicates assignment of a value to the variable. 

4. Object code listing (optional) - lists the program again with the 
generated machine language code following each statement. The 
generated code portion shows the memory type, location, hexa- 
decimal contents, and dis-assembly (assembly language equivalent) 
of each instruction generated. 

5. Subprograms called - lists name, type, and number of arguments. 

6. Statement labels - lists memory location, label name, and how 
used. 

7. Statement locations - gives source line number and memory loca- 
tion (within a block) of each non-comment and non-continuation 
line in the program. 
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8. Number of warnings and errors. 

B.2 Statement Error Diagnostics 

During compilation, statements which violate the syntactic or semantic 
rules of the language are recognized and error indications are printed. 

There are two levels of statement diagnostics: warnings and errors. 
Warnings are issued for minor infractions where the compiler can still 
determine what is to be done and compile the statement. Errors are 
severe violations of the language. In the case of errors, compilation 
proceeds as if the statement was never encountered. The statement label, 
if any, remains defined. If the error statement is ever executed, it 
will cause a link to a system routine which will terminate execution 
of the program and notify the user that an attempt has been made to 
execute an erroneous statement. The name of the program and the line 
number of the statement will be displayed. 

One character of the statement is marked with an up-arrow symbol (A) 
output directly beneath the erroneous character, for example: 



ZATA = X + Y* - A. 
A 



The character "-" is marked as an error. 

In the case of a syntax error, the marked character itself was unacceptable, 
as in the example above. In the case of a semantic error, an identifier 
or other construct is in error. The mark indicates the last character of 
the construct in error. For example, in the line: 



COMMON ALPHA, BETA, ALPHA, GAMMA 

A 



the mark indicates that the identifier ALPHA is misused 
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The compiler attempts all interpretations of statement type before dis- 
carding a statement. The marked position indicates the greatest amount 
of correct information found under the most logical assumption of 
statement type. 

A comment specifying the reason for the failure is output directly after 
the marked line. There may be more than one diagnostic per line. The 
diagnostics are listed left-to-right. Each diagnostic is followed by a 
sequence of characters: "E*E*E...E" or "W*W*W. . .W" indicating "error" 
or "warning", respectively. 

An alphabetic list of possible statement diagnostics follows: 

ARGUMENT CONVERTED (Warning) 

The type of the indicated parameter for an intrinsic function was 
converted to agree with the type required by the function. 

ARGUMENT COUNT (Warning) 

The number of parameters to a subprogram is wrong either because 
it is an intrinsic function which the compiler knows about or 
because the same subprogram was called previously with a different 
number of parameters. 

BLOCK DATA ONLY 

A DATA statement not in a BLOCK DATA subprogram attempted to 
initialize a variable in COMMON. 

An executable statement has been included in a BLOCK DATA 
subprogram. 
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CONSTANT SIZE 



The size of the indicated constant is outside the allowable 
range . 



DATA TYPE 



The type of a constant in a DATA statement does not agree with 
the type of the variables it is to initialize. 

DATA COUNT 

The number of variables in a DATA statement does not agree 
with the number of constants. 

DECLARATION CONFLICT 

An attempt has been made to declare an identifier as a FORTRAN 
entity (simple variable, array, subprogram, statement function 
name) which has already been used otherwise. 

DUPLICATE DUMMY 

A formal parameter has been declared twice in a statement func- 
tion definition, FUNCTION, or SUBROUTINE statement. 

EXTRA COMMA (Warning) 

More than one comma has been encountered at a point where a 
single comma was expected. 
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FORMAT LABEL 

The indicated statement number was declared in the label field 
of a FORMAT statement and is being used in some manner other 
than as a format reference. 

ILLEGAL DO CLOSE (Warning) 

A DO loop was closed with an illegal statement. 

ILLEGAL LABEL 

1. Statement number is less than 1 or greater than 99999. 

2. A DO statement references a previously defined label or 
a label previously referenced as a FORMAT. 

ILLEGAL NUMBER 

1. FORMAT, DATA or CONST repeat count not greater than zero. 

2. Unary minus of Hollerith or Hexadecimal constant. 

3. Illegal complex number format. 
JUMP LABEL 

1. Statement number which is not a FORMAT label has been used, 
as if it were. 

2. A FORMAT label has been previously referenced by an IF or 
GO TO statement. 



-87- 



LABEL MISSING (Warning) 

1. Indicated statement cannot be executed because it has no 
statement number. 

2. The indicated FORMAT cannot be used because it has no 
statement number. 

MISSING COMMA (Warning) 

A comma was missing at a point where one was expected but com- 
pilation could continue. 

MISUSED NAME 

An identifier has been used in the wrong contex, such as: 

- A formal parameter in a DATA or EQUIVALENCE statement. 

- A variable dimension which is not a simple formal 
parameter . 

- A subprogram name used without parameters in an ex- 
pression. 

MULTI DEFINED 

A statement number is defined more than once. 

NOT ARRAY 

An identifier which is not an array name has been used where 
an array name should have appeared. 
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NOT INTEGER 

A variable or expression of type other than integer has been 
used where only integer type is allowed. 

NUMBER OF SUBSCRIPTS 

The number of subscripts in an array reference is, incorrect. 



RANGE 



1. The second character in the declaration of an IMPLICIT range 
does not alphabetically follow the first character. 

2. A constant subscript array reference has a subscript which 
falls outside the size of the array. 

STATEMENT NOT ALLOWED 

1. A statement has been used in an illegal context. 

2. An illegal logical IF secondary statement or the statement 
is in the wrong order, such as a statement function defini- 
tion not preceding executable statements. 

SYNTAX 

1. Usually erroneous punctuation or an illegally constructed 
expression, 

2. The character marked shows how much of the statement was 
scanned before it ceased to make sense. 
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TYPE CONFLICT 



1. The same first character has been declared two different 
types in an IMPLICIT statement. 

2. The types of the operands of an arithmetic or logical 
operator are illegal. 

3. The types of the right- and left-hand sides of an assign- 
ment are improper. 



UNDIMENSIONED 

A simple variable is followed by a left parenthesis. 
UNRECOGNIZABLE 

The entire statement was unrecognizable. 
UNSUCCESSFUL COPY 

A COPY statment could not be performed. 

B.3 Program Error Diagnostics 

After the source program has been listed, summary error messages pertaining 
to the program as a whole are listed. 

The following describes each of these messages or set of messages. 

FUNCTION NAME NOT REFERENCED 

This message appears at the end of any FUNCTION subprogram in 
which the function name does not appear on the left-hand side 
of an assignment statement. 
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OPEN DO LOOPS 

Following this heading, all DO loops which were not closed are 
listed in the form: 

"statement-number OPENED AT LINE line-number" 

UNDEFINED LABELS 



All undefined statement numbers are listed after this heading. 
Each undefined statement number appears as: 

statement-number FIRST REFERENCED AT LINE line-number 

ALLOCATION ERRORS 

This heading is followed by a list of identifiers that were 

incorrectly assigned memory locations by the program. These 

errors are caused by COMMON and/or EQUIVALENCE statements. 
Such errors as: 

- EQUIVALENCE statement references variables for different 
COMMON blocks. 

- Extending a COMMON block backward. 

- Specifying an impossible equivalence group. 
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APPENDIX C 
. COMPILER OPTIONS 

C.l Introduction 

The following compiler options are placed in parentheses on the compiler 
invocation control card. 

Conditional Compilation - CO 

This option causes lines containing a D in column one to be compiled . 
Statements with a D in column one are treated as comments if this option 
is not selected. 

Cross Reference - XR 

This option causes the compiler to print an alphabetic listing of each 
identifier and where it was defined, altered, or referenced. 

List Object - LO 

This option causes the compiler to print pseudo-assembly language state- 
ments corresponding to the object code generated. 
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APPENDIX D 
FORTAN RUNTIME LIBRARY 

D.l Introduction 

The programs involved with the process of executing a FORTRAN program may 
be divided into two categories: 

Funtion Library 
Arithmetic Larbrary 

D.2 Function Library 

The Function Library consists of the basic t-x carnal functions and non- 
inline intrinsic functions. These are all referenced directly by the 
FORTRAN program, therefore, the module name is the same as the FORTRAN 
function name. 

The Function Library contains the following modules: 



ALOG 


CLOG 


DCOS 


EXP 


ALOG10 


COS 


DEXP 


IDINT 


ATAN 


COSH 


DLOG 


SIN 


ATAN2 


CSIN 


DLOG10 


SINH 


CABS 


CSQRT 


DMOD 


SQRT 


CCOS 


DATAN 


DSIN 


TANH 


CEXP 


DATAN2 


DSQRT 




Arithmetic 


Library 







The Arithmetic Library consists of the functions indirectly referenced by 
the FORTRAN program. These modules all have the prefix F% to identify 
them as part of the FORTRAN run-time library and to limit conflicts with 
user names. 
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Name 



Function 



F%CVDI 
F%CVID 
F%DADD 
F%DSUB 
F%DMUL 
F%DDIV 
F%CMUL 
F%CDIV 
F%DTNE 
F%DTEQ 
F%DTLE 
F%DTLT 
F%DTGE 
F%DTGT 
F%CTNE 
F%CTEQ 
F%PWII 
F%PWIR 
F%PWID 
F%PWIC 
F%PWRI 
F%PWRR 
F%PWRD 
F%PWRC 
F%PWDI 
F%PWDR 
F%PWDD 
F%PWDC 
F%PWCI 
F%PWCR 
F%PWCD 
F%PWCC 



Convert double precision to integer. 

Convert integer to double precision 

Double precision add 

Double precision subtract 

Double precision multiply 

Double precision divide 

Complex multiply 

Complex divide 

Double precision test .NE. 

Double precision Test .EQ. 

Double precision test .LE. 

Double precision test .LT. 

Double precision test .GE. 

Double precision test .GT. 

Complex test .NE. 

Complex test .EQ. 

1**1 power routine 

I**R power routine 

I**D power routine 

I**C power routine 

R**I power routine 

R**R power routine 

R**D power routine 

R**C power routine 

D**I power routine 

D**R power routine 

D**D power routine 

D**C power routine 

C**I power routine 

C**R power routine 

C**D power routine 

C**C power routine 
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